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We investigate a model of chaperone-assisted polymer translocation through a nanopore in a 
membrane. Translocation is driven by irreversible random sequential absorption of chaperone pro- 
teins that bind to the polymer on one side of the membrane. The proteins are larger than the pore 
and hence the backward motion of the polymer is inhibited. This mechanism rectifies Brownian 
fluctuations and results in an effective force that drags the polymer in a preferred direction. The 
translocated polymer undergoes an effective biased random walk and we compute the correspond- 
ing diffusion constant. Our methods allow us to determine the large deviation function which, in 
addition to velocity and diffusion constant, contains the entire statistics of the translocated length. 
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Proteins, nucleic acids, and various products synthesized inside cells are transported within the cytoplasm by 

^^ molecular machines or motors p] . Polymeric chains may also translocate in order to get in or out of organelles within 

the cell, or to cross the cell's outer membrane [5]. Both transport and translocation cannot rely on diffusion alone: 

macromolecules are subject to thermal fluctuations which are isotropic and from which, according to the second law 

of thermodynamics, no useful work can be extracted. Pure diffusion can account neither for the directionality of 

^I^ motion nor for the observed time scales — supplementary mechanisms that induce active and directed motion must 

be taken into account. The 'Brownian Ratchet' [3 provides a general setting to describe rectification of diffusive 

'~~^ motion using chemical energy. This paradigm helps to understand the physics of molecular motors and to construct 

'~^ mathematical models of the translocation process. As explained e.g. in Refs. [11[5], Brownian motion will cause a 

qi protein undergoing the translocation through the pore to fluctuate back and forth, and chemical asymmetries will 

rectify its displacement. Several mechanisms [5] can play this role including, for instance, disulfide bond formation, 

attachment of sugar bonds, binding of chaperone proteins on one side of the membrane. These rectification processes 

induce an effective force on the polymeric chain that drags it in the desired direction by inhibiting backwards motion. 

The chaperone-assisted translocation [U [51 [71-fT2] is a common mechanism for protein translocation and perhaps 

for DNA transport through membranes [13] . Earlier work 31 ISl 13 [13 has focused on continuous space and time 

"jrt descriptions. This assumption effectively means that the chaperone molecule is much longer than the lattice size. In 

Cj the case of nucleic acids, for instance, the lattice size is the interbase distance ~ 0.36 nm which is about 6 times shorter 

^ than the chaperone molecule. Hence the discreteness of the monomers which are the binding sites for the chaperones 

'^ can have sizable effect on the kinetics of the process. D'Orsogna, Chou, and Antal [12] have recently investigated 

^ a discrete model of the polymer translocation driven by the irreversible random sequential absorption of chaperone 

O proteins on one side of the pore. The chief analytical result of Ref. \12\ is the computation of the translocation velocity. 

I ^1 Needless to say, the translocation velocity (here, the average translocated length per unit time) provides a partial 

description which is not sufficient in many biologically relevant cases, especially when the translocated polymers are 

^~~^ not too long. In these situations dispersion (quantified by diffusion coefficients) and more generally large deviations 

^ play an important role. The goal of the present work is to carry out a thorough mathematical analysis of this polymer 

translocation model. We shall focus on the diffusion coefficient, although our approach gives the full statistics of the 

translocated length, for instance all its higher cumulants. 

fy-\ The outline of this work is as follows. In section |l| we describe the translocation model and discuss some of its 

general scaling properties. In section |ll| we investigate the simplest case when chaperone molecules attach to single 

monomers along the polymer chain. First, we calculate the velocity as a function of the parameters of the model. Then 

^^ we determine the large deviation function which encodes the full statistics of the translocated length; in particular, the 

large deviation function contains the velocity and the diffusion constant. In section III we study a continuous space 
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model where the chaperones are free to bind anywhere on the translocated segment of the polymer. This corresponds 

to the limit when the size of the chaperone molecule is large with respect to the monomer size (providing the natural 

S^ minimal size in the problem, the lattice spacing). In this continuum case the velocity and the fluctuations can also 

?— ( be determined analytically. We show that the predictions of the discrete and the continuum frameworks match when 

^ the chaperon attachment rate A is vanishingly small. In the complimentary limit of large attachment rate, the system 

becomes equivalent to a continuous time random walk. This limit admits a separate more elementary analysis and 



it provides a useful check of the consistency of our calculations. Finally in section IV we discuss possible extensions. 
A generalization to the situation when the polymer undergoes a bias diffusion and various technical calculations are 
relegated to Appendices. 



I. GENERAL DESCRIPTION OF THE SYSTEM 

Consider a polymer chain that passes through a pore in a membrane (see figure IJ). We assume that the motion 
of the polymer chain is equivalent to unbiased random walk, i.e., it is fully described by one parameter, the diffusion 
coefficient D. The pore is located at cc = and we focus on the polymer segment to the right of the membrane (the 
region located on the right-hand side of the pore is the target region of the translocation process). At a given moment, 
this polymer segment consists of a certain number L of monomer units, each of size a, labeled 1, . . . , L. Translocation 
is reversible and the polymer can move by one-monomer unit to the right or left with equal rates, resulting in a zero 
average translocation velocity and a non- vanishing diffusion coefficient D for the 'bare' polymer. We now suppose that 
the medium on the right side of the membrane contains a fixed density of special molecules — "chaperones" — that 
adsorb irreversibly, with a certain rate A, onto unoccupied adjacent sites of the polymer. A chaperone is sufficiently 
large that it cannot pass through the pore (Fig. [l}. As a result, the chaperones rectify the polymer diffusion so that 
it passes through the pore at a non-zero speed v and has a certain effective diffusion coefficient V. More precisely, L, 
the number of translocated monomers at time t, is a random variable. Let P(L, t) be the corresponding probability 
distribution. One anticipates that this distribution is asymptotically Gaussian, 

where the velocity and the diffusion coefficient of the tip characterize, respectively, the average and the variance of L 

(L) = vt, (L^) - (i)2 = 2Vt . (2) 

The quantities v and V depend on A, the length i of the chaperone molecule, the monomer length a, and the bare 
diffusion coefficient D. On dimensional grounds, wc can write 

'^ = f^(A,!), T^-DG{X,i) with A=^. (3) 

These simple relations predict scaling behaviors in the interesting regime of small adsorption rate, that is, when A ^ 1. 
In this limit we anticipate that the length of the chaperone molecules should become irrelevant. This implies that 
F ~ A^/"^ as A — )■ 0, while G remains finite. Therefore 

V - A^/^D^/^ and Vr^D when A = — -> . (4) 

In the opposite limit of large adsorption rate, A ^ 1, we anticipate that adsorption becomes irrelevant. This suggests 
that 

D - Xfi 
V '^ — and V ^ D when A — !■ cx) . (5) 



m 

(a) (b) 



(c) 



FIG. 1: Illustration of chaperon-assisted translocation, (a) The polymer can hop in either direction, (b) The polymer can hop 
only to the right because an adsorbed chaperone (large dot) is next to the pore (hole in the membrane barrier) and is too large 
to enter, (c) Adsorption of a new chaperone (shaded) at a site on the leftmost chaperone-free segment. 

In the following sections we derive exact expressions for the functions F{X) and G(A) defined in Eq. ([s]). We will 
confirm the scaling behaviors of Q-lls]) and determine analytical expressions for the prefactors. 

II. ANALYSIS OF THE DISCRETE SPACE MODEL 

In this section we analyze the simplest discrete model corresponding to the situation when chaperone molecules 
have the same length as the monomer units of the polymer chain {£ = a). Bound molecules do not overlap, that 



is, there is at most one chaperone per site. Another assumption is that the binding of chaperone molecules to the 
polymer is irreversible. No other assumptions will be made and the following analysis will involve no approximations. 

The polymer hops to the left and right with equal rates which we set equal to unity; we further set ^ = 1. With 
these agreements, the intrinsic diffusion constant D of the bare polymer is also equal to 1. Chaperones attach to 
empty sites with rate A and we assume that attachment is irreversible (the detachment rate is 0). Since we set i = 1 
and D = 1, the adsorption rate A becomes dimensionless, A = A. 

At any given moment, the complete description is provided by the length L of the translocated segment of the 
polymer and the positions of the occupied monomers by adsorbed chaperones: 

1 < mi < 7712 < . • . < rnnj < L . (6) 

The length L of the translocated segment, the number M of the attached chaperone molecules (equivalently, the 
number of occupied monomers) and their positions mi, . . . ,mM along the polymer chain are all random variables. 
At a given time, the first monomer to the right of the pore is always labeled as 1. The above definition (|6| shows 
that an occupied monomer is not allowed to move to the position a; = 0; in other words, if mi = 1 the segment can 
hop only to the right, while for mi > 1 both hops are equally possible. The full description requires the computation 
of the probability density P{L; mi, . . . , niM', t). A shorter description ignores all monomers except the one closest to 
the pore, located at a certain site j (i.e. mi = j). In other words, this shorter description focuses on the segment 
probability Qj(t), namely the probability that at time i, the leftmost chaperone is located at distance j from the pore. 
Schematically, 

Qj{t) =Prob{|o...o«} 

where the vertical line represents the position of the membrane pore. The key feature of this problem is that the 
time-evolution equations of the Qj{tys form a closed set and can be solved. (For a review of techniques developed in 
studies of adsorption kinetics, see [Ml [15].) These equations read 



dt 

dQijt) 
dt 



= Q,-i(t) + g,+i(i) + A^Qfe(f)-[2 + A(j-l)]0,(t) for j > 1, (7) 



fe>j 



Q2(t) + A^Qfc(t)-gi(i). (8) 



fc>i 



The terms in Eqs. Q are self-explanatory, e.g. the first two terms on the right-hand side describe the gain due to 
hopping, while the third gain term represents adsorption events. Equation (Is]) does not fully fit into the pattern, but 
if we add an extra variable Qoit) with Qo{t) = Qi{t), then equation ((sl) does take the same form as Eqs. Q. As 
a useful check of self-consistency one can add Eqs. ([7|-((8]) and verify that the sum J2k>aQkit) is conserved. (By 
normalization, J2k>o Qk{t) = 1 must hold for all values of t.) 

It is useful to define the cumulative variable Em{t) = X]fc>m'9fe(*)- Note that Qm{t) — E,n-i{t) — Em{t). The 
quantities Emit) represent empty interval probabilities 

^^(t)=Prob{|o_^}. 

tn 

In particular, we have E^it) = 'y.k^iiiQki't) = 1 and E_i{t) = Qo{t) + -Eo(i) = Qo{t) + 1- The equations of motion of 
Era{t) are obtained from Eqs. ([7])-([8]) by summing these equations for k > m: 

^%^ = Em-i{t) + ^™+i it) - (2 + Am) Emit) . (9) 

at 

Taking into account the boundary values, we see that equation Q is valid for all m > 0. The interpretation of the 
evolution equation ^ is simple: the gain terms arise from diffusion; the loss term describes diffusion and attachment 
(the latter occurs with rate A onto any of the m sites of the empty interval). For a more detailed derivation of 
equations 0-(|9]) see Ref. [12] and chapter 7 of Ref. [15] . 

A. Steady state weights and velocity 

In the stationary state, the values of the E„i^s must satisfy the following recursion: 

Em-l + Em+1 = (2 + Am) Em . (10) 



We emphasize again that the boundary case to = also fits into the pattern and therefore we have Ei + E-i — 2. 
Equation ( 10 1 is essentially the difference analog of an Airy equation. The solution of this difference equation that 



satisfies the boundary condition Eq = 1 can be expressed in terms of Bessel functions. This can be done by comparing 
( 10 1 with the well-known identity [T5] 



X 



for the Bessel functions. We notice that Eq. ( [lO| is solved for all values of to. (including m ~ 0) by choosing 

E„ 



:^^^^±4^ with A=| 
Ja(A) a 



Having computed these steady state probabilities, we derive the velocity of the polymer 

Ja+i(A) 



v = Qi = l^Ei 



1 



Ja(A) 



(11) 



(12) 



(13) 



This result agrees with that found in Ref. 12 . The expression ( 13 ) is compact but fairly complicated as the adsorption 
rate also enters into the index of the Bessel function. The plot oi v = v{\), see figure l2J shows that in accord with 
intuition, the velocity is a monotonously increasing function of the adsorption rate. 

To better appreciate the behavior of the velocity, let us consider the limits of small and large adsorption rate. The 
small A behavior is derived from an asymptotic formula |17j 



Jx{X + x) = 



9 \ 1/3 / r, X 2/3 

'^ Ai(0) ^' 



which is valid for A — >■ oo and for a fixed value of x. We thus obtain 

„ = -^LM>1/3 



Ai'(0)a; 



Ai(0) 



Using the well-known expressions |16j 

Ai(0) 



1 



32/3r(2/3) ' ^' ^^^ 3i/3r(l/3) 



equation ( 15 ) can be rewritten as 

V = A>}l'^ with 



-^^ ^m/sf^ =0.72901113. 



(14) 



(15) 



(16) 



(17) 



The dependence of w on A agrees with the scaling argument given in Eq. Q . This behavior can be understood more 
physically as follows. Let r be the typical time between two successive adsorptions on the leftmost empty segment. 
The tip will move a distance d ~ i/t. Now we note that Adr ^ 1 since A is the adsorption rate per site. These 
relations imply that t ~ \~'^l^ and d ~ \~^l^ . Therefore 



Ai/3 



explaining the A^/"^ scaling. Higher orders of the expansion of the velocity with respect to the attachment rate A can 
be derived from the asymptotic expansion 



Jx{X + x) 



Stt 



/ . -Sm( 



m>0 



TO,+ 1 \^/TO,+ l 

X) sm I — - — TT r 



A + x 



(18) 



extending the expansion of Eq. (14). (A proof of (18) can be found in the treatise [T7] which also contains explicit 
expressions for Bm{x), polynomials in x of order m, for small orders [18].) From Eq. (18) we deduce that 



y = AX'/''-^ + ^X'/' + 0{X') 



10 ' 140 •■ ' " V / ' (19) 

with A defined in Ec^. (17). This expression is more precise than (17l and it can be used for data fitting (see figurel2| 





FIG. 2: Velocity and difTusion constant versus the attachment rate. The left figure displays the exact formula (131 for the 
velocity, the dashed line represents the asymptotic formula (19 1. The right figure shows the diffusion constant (31 1; the small 



A expansion (dashed line) given by Eq. (|32[) is also plotted for comparison 



B. Diffusion constant and higher cumulants 



The calculation of the velocity described in subsection II A| relies on the steady-state probability distributions. The 
diffusion constant, however, cannot be determined from the steady-state distributions. Indeed, dispersion embodies 
information about transient states — it can be expressed in terms of integrals of two-time correlation functions 
|19j . while the knowledge of the steady-state characteristics is insufficient. Therefore in many interacting stochastic 
processes the analytical determination of the diffusion constant is quite involved [T9H2T] . In this subsection we describe 
an approach that leads to an analytical expression for the diffusion constant V. This approach will also allow us to 
compute the large deviation function which, in addition to velocity and diffusion constant, contains the entire statistics 
of the translocated length. 

Before embarking into calculations we notice that in Ref. pJJ, the diffusion constant has been probed through 
Monte-Carlo simulations. The value V « 0.388 was found numerically in the limit of vanishing deposition rate A <C 1, 
whereas in the A — > cxd limit the model becomes equivalent to the one-dimensional burnt-bridge model 22 , for which 
the dispersion has been calculated in [23j. We now show how to compute I'(A) analytically. 

As explained in Section |lj we are interested in the statistical properties of the translocated length L in the long 
time limit, see ([2]). Since the motion of the polymer is affected only by the leftmost chaperone molecule, let us 
disregard the positions of other chaperon molecules and focus on Pj{L,t), the probability that at time t the polymer 
has translocated a length L inside the cell and that the leftmost chaperone is at distance j from the pore. The 
evolution equations for Pj{L,t) read 



dPj{L) 
Jt 

dPijL) 
dt 



Pj.i{L - 1) + P,+i(L + 1) + A^ Pfe(L) - [2 + A(j - l)]P,{L) for j > 1, 



k>j 



P2{L + l) + XY,Pk{L)-Pi{L) 



(20) 
(21) 



fc>i 



These equations are just a more comprehensive version of Eqs. (|7|-(|8|. Using the sum rule X]l>7 ^i(^i^) = Qji^) 
which follows from the definitions of Pj{L, t) and Qjit) one can recover Eqs. (|7])-(l8]) from Eqs. (|20|-(21 ). 
We now introduce the generating function (which is essentially the discrete Laplace transform)^ 



n,(M,i)=E^''^^:'(^'*) 



(22) 



In equations (20)-(21| and hereinafter we usually suppress the dependence on time, e.g. we shortly write Pj{L) 
instead of Pj(L, t). In addition to suppressing time dependence, we shall often suppress the dependence on the 
fugacity parameter ^; for instance, IIj means T\j{pL,t). 



From Eqs. (20)-(21| we deduce the governing equations for IIj's: 



dt 

(Ml 

~df 



- = e^n,-i 



n,+i + A ^ n, - [2 + A(j ~ 1)] n, for j > i, 



k>j 



^na + A^Hfe-ni, 



(23) 
(24) 



fc>i 



Note that these equations are the same as Eqs. (l7| and (IS]) satisfied by the Qj's, except for the exp(±/i) factors. The 
infinite system of equations (23)-(24) can be rewritten in a matrix form 



dU{fi) 
dt 



Mi^l)u{^^) 



(25) 



Here n(/i) is a column vector n(/x) = (111,112, . . ■)'^ and M(/i) is a matrix 



/-I e^^ + A 



M(ai) 



A 



^(2 + A) e-^ + A 



A 
A 







-(2 + 2A) e-^ + A A 



\ 







V 



(26) 



For II — 0, the matrix 
non-vanishing /x, we can write 



e^ -(2 + 3A) e~^ + A A 
is the continuous-time Markov matrix that governs the evolution of the Qj's. For 



/ 



oMi' 



= ^nj(^,t)-e^(^)* for t^oo 



(27) 



with U{fj,) being the largest eigenvalue of matrix M(/i). (The existence and the uniqueness of the maximal eigenvalue 
is guaranteed by the Perron- Frobenius theorem [H].) In the long time-limit, the function t/(/J,) generates all the 
cumulants of L. Indeed, taking the logarithm of Eq. (27), and expanding it with respect to /i, one obtains 



.,.,^,.<^.,.=<^^y^ 



= jiv + ijl V 



(28) 



Hence, C/(/x) contains the full information for the long-time statistics of the translocation length; in particular, the 
second-order term in the Taylor expansion provides the diffusion constant. 

An implicit formula for C/(/i) can be derived following a path very similar to that used in the previous section to 
calculate the velocity v. The computations are more cumbersome [see Appendix K] but straightforward. We obtain 



U{p) - (e^ - 1) 



A _ ■^A[1+C/(m)/2] + i(A) 
V Jm1+U{u)/2\{^) 



with A 



A 



(29) 



A perturbative expansion for small ^i allows us to calculate the cumulants of the translocated length L. At the first 
order we have: 



f/(/i) = A* 1 - 



Ja+i(A) 
Ja(A) 



o(m^ 



(30) 



which agrees with the already-known expression ( 13 1 for the translocation velocity. Developing to the second order 



allows us to derive a formula for the diffusion constant of the polymer chain 

V V ( d Ja+i{x)' 



V 



X\dA Ja{x) 



x=A 



(31) 



We emphasize that the derivative is with respect to the order A of the Bessel function and not with respect to its 
argument x, i.e. the value x = A is to be substituted after the derivative has been taken. Using the asymptotic 
expansion (18 1, we deduce 



0\J 



(32) 



with A defined by Eq. (17 1. The exact formula (31 1 and the small A expansion are represented in figure l2J In the 
leading order, we find 

r>~yl^ ^3( ^f''t\ ] =0.3874382381... (33) 

The Monte-Carlo simulation results fT^ are in excellent agreement with this analytical prediction. At first sight, the 
agreement even seems too good: Equation (32) shows that the corrections to the limiting value are of the order A^^"^, 



so they could be significant even for small values of the attachment rate A. However, the amplitude in front of the 
correction term is very small (numerically 0.00851). This explains why in numerical simulations |12j the dispersion 
looks almost constant for A < 1. 

Thus we have obtained the complete statistics of the translocated length of the polymer in the situation when the 



chaperone molecules have the same length as the monomer units. Our chief result, Eq. (29), is an implicit highly 
transcendental relation for /7(/x), the largest eigenvalue of matrix M(/i). The generating function C/(^) is related by 
Legendre Transform to the large deviation function of the translocated length L (see e.g. Ref. [55] on large deviation 
techniques). We emphasize again that this generating function contains all cumulants: expanding this function, 
^(m) — J2n>i ^"^' ^^ extract v — Ui, V = 2[/2, see (30), and subsequent cumulants 



i— >oo t 

U, = hm (^V4(X)(LV3(Ly + 6(L)^(L^) ^^^^ 

i— >-oo t 

Finally we note that all results of this section can be generalized to the situation when the bare polymer undergoes 
a biased diffusion; an outline of the computations and the final formulas for the velocity and the diffusion constant 
are presented in Appendix [B] 

III. THE CONTINUOUS SPACE MODEL 

A class of proteins responsible for the translocation mechanism across the endoplasmic reticulum is known as the 
Hsp-70 family. Typical chaperones belonging to that family are of 2nm in size, representing roughly 6 or 7 lattice 
sites of the translocating polymer. The size £ of the chaperones can have an important quantitative effect on the 
translocation process, because it modifies the velocity v and the diffusion constant V, as explained in Eq. ^. For 
chaperones of length i greater than the monomer size a, the discrete model of the previous section can in principle 
be solved. The resulting expressions for the velocity are not explicit, however, because a linear system of £ equations 
must be solved in order to adjust the boundary conditions |12j . Here we examine the limiting case when the length of 
the chaperone molecule vastly exceeds the monomer size. Keeping the chaperone of size £ finite we therefore consider 
the a — >■ limit when the discrete lattice turns into a continuous line. 

We set again the length £ of the chaperone protein and the bare diffusion constant D to unity. With this agreement, 
the dimensionless attachment rate \£^ / D is equal to A. Thus mathematically the chaperon molecules are segments of 
unit length that absorb onto a continuous line with rate A. Adsorption problems where the substrate is the continuous 
line are usually solvable if the analogous problems where the substrate is the one-dimensional lattice are tractable 

m- 

Let Q{x, t)dx be the probability that the left edge of the first absorbed chaperone on the right of the pore is located 
between x and x -\- dx, see figure^ Further, denote by E{x,t) the probability that the first chaperone molecule is 
more than a distance x away. The relation between these two probabilities is 



Eix,t)= f 

J X 



dyQ{y,t) (34) 



The probability densities Q{x,t) satisfy 

dQ^dHl 
dt dx"^ 
dQ ^ d^Q 
dt ~ dx^ 



/•OO 

-A(x-1)Q + a/ dyQ{y), x>l (35a) 

Jx+l 

/•OO 

+ X dyQiy), x<l (35b) 

Jx+l 



These equations can be derived directly by studying how Q{x, t) varies between time t and t + dt. One can also obtain 
them as continuous limits of the discrete equations ([7 1 and Q . 

The governing equations for the empty interval probabilities E{x,t) are 



dE _ d^E 
dt dx'^ 
dE _ d^E 
dt dx'^ 



j-x + l 

\{x- l)E-\ j dyE{y), x > 1 



A 



dyE{y), x<l 



(36a) 
(36b) 



In the large time limit, Q{x,t) and E{x,t) become stationary. Here we focus on the stationary empty interval 
probabilities E{x) = E{x,t = cxo) that satisfy 



d^E 
dx'^ 
d^E 
dx'^ 



X{x - l)E{x) + A / dy E{y), x > 1 

J X 

rx+1 

/ dyE{y), x<l 



X 



(37a) 
(37b) 



This set of integro-differential equations can be viewed as a generalization of the classical Airy equation [26] . Essen- 
tially, we need to solve (37a); the solution of (37b I will then be found by a simple integration as the right-hand side 



is known. Equation (37a) is non-local as can be seen by differentiating it with respect to x: 

d^E{x) 



dx^ 



X{x-2)^^ + XE{x + l) 
dx 



In the following, we derive an exact solution of Eqs. (37a|-(37b) 
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FIG. 3: Illustration of chaperon-assisted translocation for the continuous model. 



A. Stationary solution and velocity 

For difference-differential equations, it is difficult even to prescribe the boundary conditions; formally one needs to 
have a boundary condition on the interval of length 1, say E{x) on interval 1 < x < 2, and then one can construct 
the solution for all a; > I. For almost any choice of the 'boundary values' the resulting solution will be unacceptable 
(diverging at a; — >■ cxd rather than vanishing, not decreasing monotonously with x, etc.). Rather than proceeding 
in a formal manner, we shall use here the Laplace method (see e.g. [27j or [28] ) which is a very powerful approach 
allowing one to deal with linear differential equations whenever coefficients are linear functions of x. This constraint 
is met in our case, so we employ the Laplace method. This method tells one to seek a solution in the form 



Eix) 



due'^='Z{u) 



(38) 



where the integral is taken over a contour C in a complex plane to be specified at a later stage. 



We determine Z{u) by inserting (38) into the governing equation (37a I. The left-hand side is 
The non-local term on the right-hand side becomes 

PX+l !■ 

/ dyE{y) ^ j due 
Jx Jc 

To simplify xE{x) we use integration by parts 



e" - 1 
"^ Z{u) 



xE{x) = f due^''xZ{u) 
Jc 



ic 



due" 



.dZ 
du 



.dZ 
du 



where the validity of the last result relies on the assumption that function e^^ Z{ii) has the same value on both ends 
of the contour C. Combining above results we find that (38) is a solution to Eq. (|37a) if Z{u) obeys 



u^Z = -\Z -\-- + \ Z 

du u 

Solving this differential equation we find 

/ u"" r , e" -I 

Z = const • exp —u r + / dv 

V 3A 7o V 

If we take the contour C that goes from infinity to the origin along the ray with arg(w) 

to infinity along the ray with arg(w) = f + e (see figure l4|), then e^^Z{u) vanishes on both ends of the contour C (as 
long as < e < 7r/6). In the e — >■ -|-0 limit, the contour C coincides with the imaginary axis. Writing u — iw we 
obtain the desired integral representation 



3-iT 



(39) 



e and then goes back 



/•OO 

E{x)^A dwC{w)C{w,x) 
Jo 



(40) 



with 



C{w) 



exp ( / dv 



cos V — 1 



V 



C{w,x) — cos[w{x ~ I) + W] 

w f^ sin V 

W = W{w) = -^ + / dv-^^ 

3A ./n V 



(41) 



size {x < 1) behavior to 



and yet undetermined constant A. Equation (40) is valid for a; > 1 and can be used to determine the behavior of 
E{x) at large x (see Appendix CI). 



The above solution (p}8|-(|39p is valid for a; > 1. Using these results we reduce equation (37b) governing the small 

sin('(i;a; + W) — sm{W) 



— - A4 '"" 
dx^ 



dwC{w) 



Integrating (42) subject to E{Q) = 1 we obtain 

E = \-Bx-\A 
with $ 



dwC{w)'^{w, x) 



sni{wx + W) — sin W — xw cos W 2 

7^ ~l~ -^ 



sinM^ 
2w 



(42) 

(43) 
(44) 
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FIG. 4: Contour in the complex plane which has been used in the implementation of Laplace's method in Eq. (38 1. The shaded 



domain (with opening angle equal to 7r/6) shows the region in which the contour C can exist and leads to a well defined function 

Z{u). 



The constants A and B are determined using the fact that the probability E[x) is continuous, and therefore the two 
expressions, (40) and (43), must match at x = 1. This gives 



with 



l-B 
A 



dwWi{w)C{w) 



cos M^ + A 



im{w + W)-{l- ^)s\nW -wcosW 



We should further equate the derivatives of E{x) at the matching point. This yields 



B 

1 

with Wt 



dwyV2{w)C(w) 



w sin W — X 



cos(u; + W) — cos W -\- w sin W 



(45) 
(46) 

(47) 
(48) 



Equation (40) and Eqs. (43)-(48) determine the probabilities E{x) in the stationary regime. The knowledge of 

-i?'(0), we deduce from 



Dabilities allows us to compute the velocity. Using v 



these stationary empty interval pro : 

( [43| that V = B. [The relation v — —E'{0) is the continuous version of (13).] Equations (45)-(47) in conjunction 

with V ^ B lead to the following expression for the velocity: 



V = 



j;^dwW2iw)Ciw) 
J^ dw [Wi{w) + W2{w)] C{w) 



(49) 



This formula can be used to plot w as a function of the attachment rate A, see figure [5] and to extract asymptotic 
expansions. In the case A — )■ 0, we make the change of variable w = X^^^uj in E q. (|49l ). Taking the A ^- limit and 
keeping w finite, we find that the dominant behaviors of the expressions in Eqs. (41), (|46|) and (48) are given by 



C= 1, W 



Wi 



3(^) and W2 = Ai/3a;sin(^) 



Substituting in Eq. ( 49 1 , we obtain the leading behavior of the speed 



A^/ 



i/Jo°°du;u;sin{^) _ 3V3r(2/3) , 



/p du! cos( 



r(i/3) 



Ai/3. 



(50) 



(51) 



This result agrees with the prediction for the discrete model, Eq. (17). This agreement is not surprising — in the limit 



of infinitely small attachment rate, the discreteness of the polymer does not play any role and the system behaves as 
if it were effectively continuous. 
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The A — ^ CX3 limit is more difficult to implement (see Appendix O), yet the final result is remarkably simple 

V ^)- 2 when A -> oo . (52) 



We can understand (52 1 as follows. In the A — > oo limit, a chaperone is attached as soon as the translocated segment 



length reaches 1. The model becomes equivalent to a continuous time random walker on a discrete lattice, that can 
only make forward steps of unit length. The randomness in the problem arises from the distribution ip{T) of the 
waiting time T between two successive jumps. The distribution '0(r) is identical to that of the first passage time at 
a; = 1 of a Brownian motion starting at a; = with a reflecting boundary at x = and an absorbing boundary at 
X = 1. The Laplace transform ^{p) = (e^-^) is obtained by utilizing standard techniques [lll[29l|30] to yield 



The velocity is then found to be 



V'(p) 



1 

(T) 



1 



cosh Y^ 



(53) 



V"(0) 



= 2. 



2.U 
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FIG. 5: Plot of the velocity v of the continuous model as a function of A^'"^ using the exact formula (491. The integrals that 



appear in the numerator and denominator of expression (491 were evaluated for various values of A ranging between 1 and 10 



(represented as dots in the figure) and a smooth interpolating curve was drawn to join these points. We observe that for large 
values of A, we have w — > 2. An empirical approximation, valid for A^''^ > 10, is given by « = 2 — -^ + 7^73 ... with Ci ~ 4 and 
C2 ~ 6. 



Calculation of the diffusion constant 



The diffusion constant 2? of the continuous translocation model can be derived by utilizing the same approach as 
in the discrete case. Namely, we shall establish an implicit equation for the large deviation function U{ii). Mathe- 
matically, one has to find the largest eigenvalue of a suitable differential operato r Ci ^ which is a deformation (with 
respect to the parameter /i) of the evolution operator for E that appears in Eqs. ( |D5p . One also has to be careful of 
how the boundary conditions are deformed. 

The details of the calculations are given in Appendix [P] and the solution for U{fi) in presented in Eq. (D15). Having 
computed U{fi), we can retrieve the formula ( |49| for the velocity (by expanding U(p) to the first order with respect to 
/i). The coefficient of the /i^ term gives the diffusion constant D. The calculations arc rather involved but systematic. 
The explicit expression for V is given in Eqs. (D16) to (D17). Although the results look unwieldy, they can be used 



to derive asymptotic limits for very small or very large values of the attachment rate A. 

In the case A — >■ we make again the change of variable w — A^/'^w. Using (50|-(51) we find the leading behavior 
of the functions 6C, SWi , SyV2 (see Appendix Id] for definitions) 



dC = 



Ai/3' 



(5>Vi = 



A 

'A1/3 



uj sm 



(5^2 ^Auj'^ cos( 



(54) 
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From this information we extract the limiting behavior of 2? by analyzing each contribution in Eq. (D16|: 



2? - 1 - « ^J° , -^+ ^J° , — +ol 55 

J^ dw [Wi + W2 C J^ dw Wi + W2 C ^ ^ ^ ' 



Substituting the dominant behaviors found above, (50)-(|51| and ( |54[ ), we arrive at 

V~l+A^ ^ ^ - ^^^^35 ^ . (56) 

Jo duj cos^ Jo duj cos%- 

Calculating the integrals in the above formula (the numerator in the third term is determined using analytic contin- 



uation and we find that the third term is equal to -1), we obtain, in agreement with the result of Eq. (33) for the 
discrete case, 

V~A^ when A ^ . (57) 

The case A — )■ oo is analyzed using techniques and tricks similar to those involved for deriving the large A asymptotics 
in the discrete case (see Appendix |D] for details). A precise investigation of the asymptotic behavior of each of the 



terms that appear in the formula (D16) for the diffusion constant can be carried out. The calculations are systematic 



and in principle straightforward, but very lengthy; the final result is shockingly compact 

2 
2? — >• - when A — )■ oo . (58) 

o 

This value is in perfect agreement with the one obtained for the continuous time random walker model with waiting 



time T distributed according to the law ipiT) with Laplace transform ')p{p) given in Eq. (53). For this effective 
problem, the diffusion constant is found to be 

2 (T)3 2 [-^'(0)]3 3' 

We remark that for very large values of the attachemcnt rate A, the chaperon-assisted translocation model becomes 
equivalent to the burnt bridge model with periodic bridges |23j . The velocity and diffusion coefficient for this model 
are known and using these results it is possible to determine the velocity and the diffusion coefficient for the chaperon- 
assisted translocation in the A — > oo limit. The connection with the burnt bridge model with periodic bridges holds 
even in the general discrete case with arbitrary £/a [12]; if i/a — )■ oo, the limiting values v = 2 and V = 2/3 
are recovered. Finally we note that the generating function technique employed in this study can be adapted to a 
calculation of the large deviation function of the burnt bridge model. 

IV. DISCUSSION 

The mechanism that drives the translocation of a protein through a membrane depends on chemical asymmetries 
between the cis and trans sides of the membrane [1]. There can be several mechanisms leading to asymmetry and 
they can be simultaneously present. In this work, we have studied a model where asymmetry arises due to the 
preferential binding of proteins on the translocated segment of the polymer. The mathematical effect that gives raise 
to directed motion is the Brownian ratchet model that plays a crucial role in the field of molecular motors [5TH5^ . It 
is important to note that the ratchet paradigm is quite insensitive to the precise mechanism at the molecular level. 
The key ingredients are the breaking of spatial symmetry and a supply of free energy which insures here a sufficiently 
strong binding of the chaperone on the polymer. In biological context, the free energy source is usually provided by 
ATP hydrolysis 0135]. 

From the point of view of statistical physics, it is interesting to note that these models, which are effectively 
one-dimensional, can be mapped to variants of driven diffusive gases and can be studied through techniques which 
were developed to understand the physics of non-equilibrium systems. These methods not only reveal underlying 
time-reversal symmetries of the system, but also provide powerful computational tools. For example, the technique 
based on the deformation of the master operator can be traced back to the proof of the Gallavotti-Cohen fluctuation 
theorem for Langevin dynamics and for Markov processes by Kurchan [36| and by Lebowitz and Spohn |37j . In the 
present work, we have assumed that the attachment of chaperones is irreversible; if the detachment rate were positive, 
the fluctuation theorem would be valid. The analysis carried out in this paper heavily relies on the fact that the 
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equations of motion depend only on the position of the chaperone closest to the pore. If the chaperones could detach 
from the polymer, the above property would no longer be true and the positions of all the chaperones bound at a 
given time to the polymer must be taken into account. The system would turn into an M— body problem with both 
M (the number of currently attached chaperon molecules) and L (the current length of the translocated polymer) 
being random variables. This appears to be a very challenging unsolved model. Besides, similar descriptions appear 
in mass transport models within a fungal hypha with mutually excluding particles progressing in a stochastic manner 
along a growing one-dimensional lattice |38| . 
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Appendix A: Calculation of the Cumulant Generating Function 



In this Appendix, we show that the maximal eigenvalue U{ji) of matrix M(/i), defined in Eq. (26), is indeed given 
by Eq. m. 



Let Gi, G2, . . . be the components of the eigenvector corresponding to the maximal eigenvalue U{^). The eigenvalue 
equation reads 



[/(m)Gj = e''Gj_i+e-^G,+i+A^Gfc-[2 + A(j-l)]G, for j > 1, 



k>j 



C/(/i)Gi = e-''G2 + A ^ Gfe - Gi = e^Go + e-^Gs + A ^ G^ - 2Gi ^ 



(Al) 
(A2) 



k>l 



fe>l 



where, in order to fit the equation for j = 1 into the general pattern, we have defined e'^Go = Gi. As in the beginning 
of Section im it is useful to introduce r,„ = X]fc>m^fe' defined for m > —1. Then, for all tti > 0, the following 
equation is satisfied: 



C/(/x)r™ = e-'^Tm+i - (2 + Am) T™ + e^E^^i 



(A3) 



The value of Fq is fixed by imposing the normalization constraint: Fq = 1. (This overall proportionality constant 
does not affect the final result.) Equation (A3 1 for to = then gives 

(A4) 
(A5) 



e-'^Fi + e^F_i = 2 + U{fi) 
In terms of the F's, the boundary condition e^Go = Gi becomes 

F_i+Fie"'' = 1 + 6-'' . 



To solve Eq. (A3 1 with the boundary conditions (A4) and (A5), we use the following ansatz: 



Jr(x) 



,mM 



where r and x have to be determined. Combining Eqs. (A3) and (A6), we get 



Hence 



a; = 2/A = A and rA = 2 + U{n) . 



^ ni. -.-/IN ^ 



MA) 



The normalization Fo = 1 is satisfied. We eliminate F_i between Eqs. (A4) and (A5) and obtain 



Fi(l - e-'") = p/ - 1 - t/(/i) = e^ + 1 - rA . 
Taking Fi from (A8) and inserting it into ( |A9| ) we obtain 

Jr+i(A) 



rA+(e^-l)- 



Jr(A) 



e^ + 1. 



(A6) 

(A7) 

(AS) 

(A9) 

(AlO) 



Using C/ = rA — 2 and extracting rA from (AlO) we arrive at the announced implicit equation (29 1 for t/(/i) 
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Appendix B: The case of a driven polymer 

Here we generalize the analysis of Section |ll] to the case of an asymmetric motion of the translocating polymer. We 
suppose that the polymer is driven by a non-zero force F through the pore: this can be modeled by saying that the 
polymer hops to the right with rate p and to the left with rate q with p/q = e^ . The basic equations now become 

^ = pQj^i + qQj+i + A ^ Qfe - [p + g + X{j - 1)] Q, for J > 1, (Bl) 

^ = gg2 + A^Qfc-pQi. (B2) 

k>l 

One has to define Qo such that pQo = qQi- 

In the stationary state, the equations for the i?,„'s defined as £,„ = J2k>m Qk become 

pE.,n-i + qE„i+i ^ {p + q + Am) E^ , (B3) 

with Eq = \ and pE^i + qEi = p + q. The solution of this difference equation that satisfies the boundary condition is 

E.^^(?-y^^^^±^ with A = ^ and A = ^^ A . (B4) 

\qj Ja(A) a p + q ^ ' 

This leads us to the velocity of the polymer: 

v=p + qiQ,-l)=p-qE,=p-^M^^^^±^, (B5) 

Similarly, we derive the following equation for the maximal eigenvalue and the corresponding eigenvector for the 
cumulant generating deformed matrix: 

pe'^r„_i + qe-tTrn+i ^{p + q + Xm + U{p)) r„ , (B6) 

which is valid for m > with boundary conditions Fg = 1 and p(r_i — 1) + ge^^(ri — 1) = 0. This leads us to the 
implicit equation for U(fi): 



Expanding (B7) to the first order in /i we recover velocity and an expansion to the second order gives the diffusion 
constant: 



^ -^ ^ V f d Jk+i[x) 



2 "^^^ A \dK Ja{x) 



(B8) 

x=A 



Appendix C: Asymptotic expansions for the continuous polymer model 
1. Large x behavior of E{x) 



To compute the large x asymptotic, we apply the steepest decent technique on Eqs. (38H(39). These equations 
show that we need to integrate e^^^' along the contour in the complex u plane with 

,,3 /•" p" - 1 

F{u) =u{x-l)-^+ dv (CI) 

3A Jq V 

We must find the stationary point u.^, where F'(u.^,) vanishes and deform the contour so that it passes through the 



stationary point in the direction of steepest decent. Using (CI) we get 



x-l- ^ + = (C2) 

A w* 
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for the critical point and we see that the direction of the steepest decent is along the imaginary axis. When a; ^ 1, 
we have 



X* = -y/X{x - 1) 



1 



2(.T - 1) 



0{x 



^5/2n 



(C3) 



Writing u = Us,+iw and knowing that the constant that appears in (39) is equal to (— iA/2), we recast (38)-(39) into 

A 



E^ -e 
2 



F, 



dw exp 



X- 1 
A 



(C4) 



where i^* = F{ui,). Computing the Gaussian integral in (C4) and using (C3) to simplify F^ we deduce 



E = A- 



o7-B, 



-A3/4(a;- 1)1/4 



exp 



VA(x-1)3A 



where 7e = 0.5772 ... is the Euler's constant. 



2. Large A behavior of the velocity 



We calculate the leading behavior of the velocity in the A — )• oo limit by studying separately the numerator and the 
denominator in Eq. (49). Keeping only the dominant terms in A, we find, using Eqs. (46)-(48l, 



dwyVi{w)C{w) ~ A / dw 



C(w) 



w 

oo 



sin(u; + W) 



l~ —\smW -wcosW] +0(1), (C5) 



dwW2{w)C{w) ~ -A 



dw 



C{w) 



w^ 



{cos{w + W)~coiiW + wsi\\W) +0(1). 



(C6) 



The function W{w), see (41), is given by 

W{w) 



^^ — +3a"^°(")+3A 



where we have denoted by Wo{w) the leading behavior. One could expect that the dominant behavior can be found 
by substituting Wq in place of W in the expressions (C5) and (C6). But in fact, we have 



dw 



Ciw) 



('sin(w + Tyo)- (1- ^)sinM^o-wcosl^o) = 0, 
( cos(w + T4^o) — cos W^o + wsin VFo) = 0. 



C{w) 



(C7) 
(C8) 



We now explain why these integrals vanish identically. Because the integrands are even functions of w (recall that 
C{w) is even and Wo{w) is odd), we can replace the lower bound in both the integrals by — oo. 

Before deriving ( |C7[ )-(C8), we shall prove a much simpler identity. Let us show that the following relation is 
satisfied: 



dw C{w) cos(w + Wo) — - 
2 



dw C{w) cos(w + Wo) = . 



We can add to this expression the following integral which vanishes because the integrand is odd 

dw C{w) sin(w + Wo) = . 



Then, using the explicit expressions (41) for C{w) and Wo(u;), and defining the function 



fw iv _ ^ 

£{w) = I dv 



(C9) 



(CIO) 



(Cll) 
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the identity ( C9 1 becomes equivalent to 



du;e^("')e™=0. 



(C12) 



This integral converges in the w-complex plane for w with positive imaginary value (i.e. 9(w) > 0). The value of 
the integral will therefore be unchanged if we deform the real axis (— c», +cxd) into any contour located in the upper 
half-plane. One can take for example w = y + ib with — cxd < y < oo and a fixed, strictly positive, value of b. 

~^ appears as a prefactor and that the rest of the 

A more 



Substituting this expression for w in Eq. (C12), we observe that e 

+00 limit proves that the value of integral is indeed 



integral remains bounded as b varies; taking the b 

elegant way to prove that the integral in Eq. (C12 1 vanishes identically is to deform the real axis as shown in Fig. [6] 




FIG. 6: Contour deformation in the upper complex plane used to prove Eqs. (C5l, (C7| and (C12l. The final contour follows 



the positive imaginary axis down from -|-Joo to and then back from to +ioo. There are no singularities or cuts on the 
imaginary axis, therefore the two halves of the contour simply cancel with each other. 



Utilizing the same approach, one can prove the following identities: 



dwe^^""^ 





poo 


dw 


e^ 


[w) 


gjtu 


- 1 




J — 00 








w 




OO 


dwe^^""^ 


e™ 


- 1 - 


iw 








W2 




00 












(w 


)^™ 


-1 


- 


iw 


+ wy2 



0, 
0, 
0. 



(C13) 
(C14) 
(C15) 



Note that Eq. (C8I is the real part of Eq. (C14) and Eq. (C7| is the imaginary part of Eq. (C15) 



In order to find the asymptotic expression for the velocity in the large A limit, it is thus not enough to replace W 
by Wq in the expressions (C5) and (C6) and we must perform an expansion to the next order. The leading behavior 



of the integral / dw Wi (w) C{w) is obtained by subtracting from the right-hand side of Eq. ( C5 ) the left-hand side of 
Eq. (C7|. There are three terms which wc now analyze separately. For the first term we obtain 



dw 



C{w) 



w^ 



{sm{w + W)- sm{w + Wq)) 



/•OO 

2 dw 
Jo 



A2/3 



C{w) . fw^\ 
w^ 

dx 



sm|-jcos(^^+g^ 



sm 



6 



Wo 

^3 



°° C(Ai/3x) a;3 / ,,3 x^ 
^ ' "- — cos A^Z-^x + — 



6 



Wo{\^'^x) 



where in the second line we have made the change of variables w —>■ X^'^x. For large values of the argument w, we 
have C{w) ~ e^'^ jw and Wq[w) — ^ 7r/2. Therefore, the above term is of order ©(A"^). 
For the second term, we find 



~*„^(i 



-) (sin(T4^) - sin(T4^o)) = 



„.S^(.'- 



dx 



Ai/3 



sina;'^/6 



A2/3 
2 



x^ ( x^ 
sm -- cos — - 
6 V 6 



Wn 
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where in the last equation we have inserted the asymptotics of C and Wq for large values of the argument and kept 
only the dominant term. 
Finally the third term reads 



oo {^ ( in\ 9 

'^^ ^;^ (cos Wo - cos W)^—j^ 
a w ^ ' Jo 



dx —-^ — - sm — sm — + Wo{X^^^x) 
x'^ V 6 



and it scales as A ^'"^ since C{w) ^ w ^ for large w. 

Comparing these asymptotics we see that the leading contribution arises from the second term. Thus the integral 



( C5 ) is equal to 



dwWi{w)C{w) = ~X^^^e'<- r dx (^21^^11^ 
Jo \ ^ 



-A2/3 



p7E 31/6 r(2/3) 



in the leading order. Using a similar reasoning, wc find the leading behavior of the integral (C6) 



/"OO />oo 

I dwW2{w)C{w) = 2X^/^6'"' dx 
Jo 



sina;'^/6 



(C16) 



(C17) 



Substituting (C16)-(C17) into the formula (49) for the velocity shows that v — > 2 when A — > cjo as written in Eq. (52 1 



Appendix D: Generating function of the continuous polymer 



In this appendix, we investigate the time evolution of the probability density P{L, x, t) in the continuum setting. 
By definition, P{L,x,t)dLdx represents the probability that the total length of the translocated polymer inside the 
cell is between L and L + dL and that the distance of the leftmost adsorbed chaperone from the pore is between x 
and x + dx. The probability density P{L,x,t) evolves according to 



dP d^P d^P d^P 

= 1-2 1 

dt dx'^ dxdL dL"^ 



/•oo 

-A(x-l)e(x-l)P + A / dyP{L,y,t). 

Jx+l 



(Dl) 



The Heaviside function Q{x — 1) ensures that this equation is valid for all values of a; > 0. If we study what happens 
in the vicinity of the a; = boundary, we observe that the correct boundary condition is given by 



dP dP 

dx dL 



= for a; = 0, L > . 



(D2) 



(This boundary condition can be obtained by writing the equation in a discrete form near a; = 0.) It is useful to 
consider the Laplace transform of P(i, a;) with respect to L 



P{^i,x,t) 



/•OO 

/ dLe''^P{L,x,t), 
Jo 



and to use the cumulative function 



f(/x,a;,t) 



dyP{tJ-,y,t). 



This function £ satisfies the evolution equation 



d£ _ d^£ 
dt dx'^ 



- X{x - l)Q{x - l)£ - X / dy£{^i,y)-2^l—+^JL''£ = C^,£, 

Jmax(l,a;) OX 



as it follows from (Dl|. The boundary condition is deduced from Eq. (D2| to yield 

92 £■ d£ , , ^ 



(D3) 



(D4) 



(D5) 



(D6) 



We observe that equation ( D5 ) is very similar to equations ( 36a 
the only difference being in the terms containing /i and /i^. T 



)-(36b I satisfied by the cumulative probability E{x, i), 
le differential operator C^ that governs the evolution 
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of £ is a deformation with parameter /i of the evolution operator Cq for E. The long time limit of £ will be given by 
the dominant eigenvector which corresponds to the largest eigenvalue C/(/i) of C^. We denote Ef^{x) this dominant 
eigenvector. It satisfies 



([/(m) - m') e. 



<PE,. 



^^2 A(a; - 1)9(2; - 1)^;^ - A 



x+l 



max(l,a:) 



dyE^,{y) -2^ 



dE^ 
dx 



The boundary condition ( D6 ) becomes 



E''{0) = ^iE'M, 



(D7) 



(D8) 



where the prime denotes the derivative with respect to x. We also normalize £^^^(0) — 1, which is an overall constant 
that does not affect the final result. We emphasize that this situation is perfectly analogous to the discrete case where 
the Markov operator M(0) was /i-deformed into M(/i) [see Eq. (26)] and where the dominant eigenvalue allowed us to 
determine the cumulant generating function. 

For a; > 1, we obtain, using again the Laplace method [571 [2S], an integral representation for E^{x), which assumes 
a form analogous to the one found in Eq. ( 40 ) 



/■oo 

Ef,{x) ^ A^ I dwCf_,{w) Cf_,{w,x) 
Jo 



Cf,{w) = exp 



dv 



cosu - 


-1 


A* 2 

-a" 


V 





(D9) 



Cf_i{w,x) — cos[u;(a:: — 1) + W^ 



W^{w) = W{w) + 



Ujp) /i^ 
A 



w 
3A 



, smw U(n) — n 

dv 1- ^ — —w . 

V A 



We shall write W^ instead of W^,(ti;) in equations below. Note that ioi fi = we recover the expressions in Eq. (41 1 
The constant A^ is yet to be determined. 

In the < a; < 1 region, the solution is given by 



r+x I ^r^x 



E^.{x) = 



gr+x _^ g 



-s„ 



-AA 



dwC'^{w)^f,{w,x) 



(DIO) 



with r± = r± (/i) = /i ± \JU{jj) and 



$^(w,a:) 



w — r+r_ 



2/Ll 



9 ^^ / 9 9s SmiwX + Wu) - T^ ^^T^ ^ C0?,(WX + W n) 

r-2, _L ,,,2^7.^ j_5/i2^ ^ t^' (ri + w^)(r_ + w^) 



wivj^ + 'uP'')(r_ + w^) 



(Dll) 



+ 



w(rj^ — rS) 



cos W,, 



r+(r^ + w^) r_(r^+u;2) r_ r^ 



^r+x 



rt_ + w^ 



Note again the similarity of structure between these equations and (43)~(44) 



The constants A^ and i?^ can be determined by taking into account that E^ix) and its derivative are continuous 
at a; = 1. These constraints lead to 



Ap,Ji^fi — 



Ai_iJ2.jj. — B^ 



2 '^ 



''4 



- r_ 

r_|_e'"+ + r_e'"~ 
2 



(D12) 



where we have used shorthand notation 



Jp.jj. — 



dwWp,^{w)C^{w), p=l,2 



and 



Wi,^(u;) = cos W^^(w) + A $^(w, 1) 
W2,^(w) = w smWfj,{w) - A 



(9$^(u;,x) 



dx 



x=l 
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Equations (D9|-(D12) determine the eigenfunction E^{x) in terms of the eigenvalue U{ii). For /i = 0, U{Q) — and 
it can be verified that E^{x) is identical to E{x) given by Eqs. (45|-(48). 

The maximal eigenvalue U{^) is implicitly determined by the self-consistent equation obtained by substituting 
X = into Eq. (D7), imposing the boundary condition (D8|, and using the normalization E^[Q) = 1. One gets 



U{fi)^fi^~fiE'{0). 



A'w,x) 



dx 



Using relation 

£"' (0) = /i — Bfj^. Inserting this into (D13) we simplify it to 



x=0 



[which can be deduced from (D12[)] and specializing Eq. (DIO) to 



(D13) 
we find 



c/(/i) = ^^B^ 



(D14) 



Solving (D12| and comparing with (D14) we arrive at an implicit expression for the dominant eigenvalue C/(/i), 



Uifi) = fi 



(r+e'-+-r_e'--) Ji^^ + (e'-+ 



I J; 



2./X 



which, using that r-|_ — r_ = 2-\/t/(/i), can be rewritten as 



Uifi) 



(r+e'-+ + r_e--) Ji,^ + (6^+ + e^-) J2,, 



n 2 



(r+e'-+ -r_e'--) Ji,^ + (e'-+ 



)J2,, 



(D15) 



We know from Eq. (28) that C/(/i) generates the cumulants of L in the long time limit: C/(^) = jjlv + ji^T) + . . .. 
Due to the presence of the overall factor ji in Eq. (D15), we observe that the expansion of C/(/J.) with respect to // can 
be generated order by order in a self-consistent nranncr. At lowest order, taking the limit /i = on the right-hand 
side of Eq. (D15), we retrieve the formula (491 for the velocity v. At the next order, retaining the terms linear in /i 



in r_|_, 7'_, Wi^p, W2,p and C^ (and substituting U{^i) — > /if) we obtain a formula for the diffusion constant T): 



V ^ 1 



^1^3 A2-i;(Ai + A2) 

" 3 "^ /p°° dw [Wi (u;) + W2 (w)] C{w) 



(D16) 



with 



Ai = / dw [yVi{w) 5C{w) + (5Wi(w) C{w)] 
Jo 

/"OO 

A2 = / dw [W2{w) SC{w) + dW2{w) C{w)] 
Jo 



The functions SC{w), SWi{w) and SW2{w) represent the first order derivatives of the functions C{w), 'Wi{w) and 
yV2(w) with respect to the parameter /i. For instance, SC{w) — — fii^l 



d^ 



and its explicit expression reads 



SC{w) 



w 

T 



Ciw). 



The explicit expressions for SyVi{w) and 6yV2iw) arc quite cumbersome 



SWi{w) 



A 



w sin W - 



A 

w 



■ sin(w -t- W) + s'mW 



- cos(w + W) 



Xw'^ 



cos W 



V 
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1 + 



V 

Xw 

v/6 



2 + w/2 



w^ 
V 2 



- V 



2A 



5W2 (w) ^ -w^cosW 

A w 



{Mw + W){^-^ + ^ysinw{l + l + 



V 



cos{w + W) 



sin W I 1 
cosM^ 



Xw-^ w'^ 
V V 2 
6 



)} 



Aw 



2 + w/2 



V 



)} 



W^ \A W W-^ 

These formulas can be used to establish the limiting behavior of the diffusion constant when A is very small or very 



large [see Eqs (57) and (58) 
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